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We present a set of scalar master integrals (Mis) needed for a complete treatment of massive two- 
loop corrections to Bhabha scattering in QED, including integrals with arbitrary fermionic loops. 
The status of analytical solutions for the Mis is reviewed and examples of some methods to solve 
Mis analytically are worked out in more detail. Analytical results for the pole terms in t of so far 
unknown box Mis with five internal lines are given. 



I. INTRODUCTION 



(N 

■"si" ■ Bhabha scattering is proposed to be measured at the International Linear Collider (ILC) in the very forward 

' direction with a precision that would allow to determine the luminosity with an accuracy of 10~* y,Q. The 

' Monte Carlo programs which were in use for the analysis of LEP data aimed at a slightly lower accuracy, and one 

CIh^ could neglect certain two-loop corrections there; for a discussion see At the ILC, the theoretical prediction of the 

ji^,. differential cross section for Bhabha scattering 

e+(p2) + e-(pi) ^e+(p4)+e"(p3) (1) 

^ ^ has to include the complete virtual photonic two-loop corrections. This is a highly nontrivial task, but with substantial 
, recent progress in several directions. Here, we will concentrate mainly on efforts to determine the two-loop corrections 
5^ ■ in a calculational scheme with a finite electron mass m, and with a regularization of both UV and IR divergences 
\ with dimensional regularization |78| . Because the Monte Carlo programs for the treatment of real bremsstrahlung 
assume a finite photon mass at intermediate states of the calculation 0, IS B , o ne will have to care about this 
fact when results will be finally combined. This might be done similarly as in |10lllll|. Alternatively, the soft photon 
bremsstrahlung can be recalculated completely in D dimensions as in |l2j| . where this is done for the simplest subset 
of the corrections. 

In the QED model with three leptonic flavors one has 154 Feynman diagrams (with all the IPI diagrams, but 
without loops in external lines), among them 68 double-box diagrams '13^. Besides the usual problems of efficient 
bookkeeping, the main problem is the evaluation of the loop integrals. One has to solve Feynman integrals with L < 2 
loops and N < 7 internal lines, 

afx) = ^ I d^k,...d^kL X 

where X — 1, kia, kiak2p, ■ ■ ■ stands for tensors in the loop momenta. This might be done by a procedure with three 
subsequent steps: 
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(i) reduce all tensorial loop integrals to scalar integrals, 

(ii) reduce these to a smaller set of scalar master integrals (Mis), 

(iii) evaluate the Mis. 

A completely numerical a ppr oach might also be possible 0| . For checks in the Euclidean region this has been proven 
to be a powerful tool 0,^3; see Appendix El 

Step (i) may be considered to be solved by now, the second one is solved for Bhabha scattering with this article 
|79j| . and step (iii) is solved for all self energies and vertices, but remains largely unsolved for the two-loop boxes. 
So, the evaluation of the two-loop boxes is the remaining bottleneck. Indeed, very few of the 33 double-box master 
integrals have been determined completely [itLITs. 19, 20] or to some extent 0|. 

The article is organized as follows. In Sectional the diagrams and prototypes for two- loop Bhabha scattering are 
identified and the method for their evaluation in terms of Mis is outlined. We discuss in Section IIIII the complete 
set of Mis and give an overview with figures and tables. The status of analytical (or semi-analytical) solutions of the 
known Mis is summarized. Some typical techniques for Ml evaluation are demonstrated. The section includes also a 
discussion of Mis with numerators. We close with a Summary. Appendix ^ includes an extension of the numerical 
method for the evaluation of Feynman diagrams with sequential sector decompositions, which allows to treat integrals 
with irreducible numerators. The final Appendix IbI defines a complete and compact set of Harmonic Poly logarithms 
up to weight 4. 



In a theory with only electrons and photons, one has in the 't Hooft-Feynman gauge 52 IPI two-loop Feynman 
diagrams, all of the double-box type. After adding to this number the IP-reducible diagrams without loop insertions 
in external lines, there are 94 diagrams. We will call this the one-flavor case Due to the existence of up to 

two closed fermion loops in certain diagrams, a complete picture of the process arises when besides the electron two 
additional flavors are taken into account. Then, there are 68 double boxes and 154 two-loop Feynman diagrams in 
total 'si'l. Fortunately, there are much less two-loop structures to be calculated; they are represented by prototypes. 
Prototypes are irreducible (sub)diagrams of a certain topology with account of the various, different masses of internal 
lines. All the Feynman integrals G{X) (see (j^J) with the same topology and propagators, but with arbitrary powers 
i^i > of these propagators, and potentially with irreducible numerators X correspond to one prototype. In Figures 
n to 121 we show the sets of two-loop self-energy, vertex and box diagrams for which the scalar Mis are needed. As 
mentioned, there are two bosonic self energies SEl and SE2, and the other three fermionic self energies renormalize 
the external lines. Further, there are five two- loop vertex and six two- loop box diagrams for Bhabha scattering. 

Technically, one has to calculate all the Feynman integrals G{X) related to these diagrams by a reduction of integrals 
with irreducible numerators X and denominators with higher powers Vi to a smaller set of scalar master integrals. 

We have determined such a set of scalar Mis for the virtual two-loop corrections to Bhabha scattering with the 
Laporta-Remiddi algorithm pH |23 . Realizations of the algorithm are SOLVE the Maple program AIR and 
the C++ library DiaGen/IdSolver |25|. We use DiaGen/IdSolver [s^, which allows to tackle two problems: (i) derive 



II. FROM DIAGRAMS TO PROTOTYPES AND MASTER INTEGRALS 




SEl 



SE2 




SE3 



SE4 



SE5 



FIG. 1: The two-loop self-energy diagrams for massive Bhabha scattering. SE3 to SE5 are needed for the renormalization of 
external fermion lines. 
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FIG. 2: The two-loop vertex diagrams for massive Bhabha scattering. 




FIG. 3; The two-loop box diagrams for massive Bhabha scattering. 



an appropriate set of algebraic equations with integration by parts (IBP) |2a| and Lorentz invariance (LI) identities 
|27j |83j |. and (ii) determine a list of master integrals by solving this set of equations. The procedure is heuristic. To 
be safe about the completeness of the solution, one has to solve the system of equations to high powers of numerators 
and denominators. If the coefhcients are kept as exact functions of the kinematical variables s,t,m? (and for several 
flavors of a further scale A/^), the algorithm is time and computer memory consuming. To minimize computer 
ressources usage, evaluation homomorphisms are used, i.e. the system is solved by projecting the coefficients to the 
field of rational numbers with suitably chosen values of the parameters. Let us mention that, from the point of view 
of complexity, the complete two-loop Bhabha scattering calculation with massless fermions is as complicated as the 
two-loop massive vertex case: on a 2 GHz Pentium PC with 1 GB memory, it takes minutes to be solved. Further, 
the number of Mis is moderate, e.g. 5 massless Mis compared to 22 Mis for the massive box diagram B3. 

There is some arbitrariness in the choice of masters. We prefer to present a set of MI without numerators. But 
we then have to allow for higher powers of propagators {vi > 1, dotted lines). This choice has one basic advantage: 
the Mis are independent of the momenta flowing inside loops. Of course, in the set of solutions for Mis, there are 
algebraic relations between scalar integrals with numerators and those with dotted denominators: so one always has 
a freedom of choice. For an explicit evaluation of the Mis, this is of importance; see the discussions in the subsequent 
sections. 

In Tabled we give a list of the net numbers of (two-loop) -I- (one-loop) master integrals needed for the evaluation 
of all two-loop vertex and box diagrams for the one-flavor case. The case of several flavors is separately discussed in 
Section UTTfI 



III. MASTER INTEGRALS 



Here we will describe the sets of master integrals needed for the calculation of all Feynman integrals for the diagrams 
in Figures ^ to |21 We use a nomenclature where e.g. V312m is the name of the integral for a Vertex with 3 lines, 
among them 2 massive lines; B514md is the name of the integral for a Box with 5 lines, among them ^ massive lines 
and one line with a dot. Sometimes there are several candidates for the same name, e.g. V614ml and V614m2 in Figure 
El or V411mldl and V411inld2 in Figure |H1 

Because they are lengthy, we give in this article practically no explicit expressions for the master integrals. Instead, 
for all the self energies and vertices, they may be found in [13j in form of a human readable Mathematica file. 
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Master sBhabha.m. We determined all these Mis, and several of the box masters, in terms of Harmonic Polylogarithms 
(HPLs) 28] . The file allows a determination of the expressions in form of polylogarithms. We used HPLs up to the 
weight 4 and give in Appendix El a complete basis for their systematic calculation. The complete list of HPLs may 
be found in the Mathematica file HPL4.m in Numerical checks were performed with the numerical integration 
package sectors. m, see Appendix 1X1 

We evaluated most of the Mis with the method of differential equations [^D] , which has been described in detail 
m many papers, e.g. [11 mm 113, and we will not repeat this here. Nevertheless, it might be useful to indicate 
some technical details of our calculations. The master integrals are defined in a Minkowskian metric, and the external 
momenta are introduced in With {pi +P2) = {P3 +P4) and pf = m? , we define s = [pi +P2Y, t = {pi — pa)^, u = 
{Pi ~P4)'^- The analytical results are expressed by dimensionless variables x and y, 



(3) 



V-s + 4 - 



corresponding to s = —(1 — x)'^/x, and y is obtained by replacing s by t. We set the electron mass to unity, — 1. 

In the rest of the article, we discuss the set of master integrals and their expansions in the parameter e = (4 — D)/2. 
First we treat only electrons and photons. The diagrams with additional flavors and thus with a second mass scale 
are discussed in Section ITlI I*'! 

A. One-loop Master Integrals 

There are five one-loop master integrals needed for the evaluation of the two-loop diagrams; see Figure 01 Our 
normalization of the momentum integrals is chosen such that the one-loop tadpole becomes: 



,D/2 J q^-i e \ 2 J \ 2 3 

For completeness, we should mention that a full calculation of the Bhabha scattering process ^ also includes the 
one loop corrections in the electroweak Standard Model (plus some higher order corrections). For their treatment we 
refer to S H IM IM S 111 IS H El El El . 





TlUm SE212m SE210m 




V311m B412in 
FIG. 4: The five one-loop Mis. External solid (dashed) lines describe on- (ofr-)shell momenta. 

B. Two-loop 2-point Master Integrals 

There are six two-loop 2-point Mis; see Figure [SI The MI SE310m may be exp ressed by two subsequent one-loop 
momentum i nteg rations in terms of F functions whose e-expansion is trivial |44l |45| . Explicit expressions may also 
be found in quoted from ^40]. The Mis SESllm and SE313m for the renormalization of external fermion legs 
are needed only on the mass shell; they may be calculated with 0N-SHELL2 with an appropriate change of 
normalization with respect to JH). For MI SESlSm see also [Uli^. The MI SE513m was first calculated to 0(1) in 
13 Hi (see also ^) and to 0(e) in IHIIl- The two remaining Mis SE312m = SE312m(l, 1,1,0) and SE312ind = 
SE312m(l, 1,2,0) are defined by: 



d^kid^k2 {kik2)- 
[{ki + ka-pY-m?]^ [^lY [kl 



SE312m(a,5,c,d) = =-^7^ / 777—^^ ^2ih r,.2i„ r,.2 ^2ic - (5) 
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These two Mis have been expressed in in terms of Tlllm and the functions SE312m(l,2,2,0) and SE312m(l, 1,3,0), 
and the e-expansion is determined there (after equation (15)) in terms of polylogarithms. 

An equivalent result, in terms of HPLs, is given in j^^i where the integrals SE312m and SE312mN = SE312m(l, 1,1,-1) 
have been chosen as masters. By an algebraic relation, valid for — 1 and — s, 

-(l + s) + e(2 + s) 2(1 -e), , „ , 

SE312md = — ^ i SE312in+ (Tlllm^ + 3 SE312mN) , (6) 

one may derive then SE312md. 

For a direct determination of the Mis with differential equations, we constructed a differential operator for the self 
energies using the scaling property: 

SE(Ap2, Am^) = A^™[S^(P'^"'')1 SE(p2, m^), (7) 
where dim(SE) is the dimension of the 2-point MI SE. The differential operator is: 

d d 
s— SE(s, m^) = -TO^— ^SE(s, m^) + dim[SE(s, m^)] SE(s, m^). (8) 
OS om'^ 

With this operator, one may derive coupled differential equations for the two masters and solve them with account of 
boundary conditions at the kincmatical point s = 0. Because the integral SE312md is one of our masters, we reproduce 
it here explicitely: 

SE312md(x) = -L + i__(^i_^ + i±|i/[0,:,] + 2±|li7[0,0,x])+O(e). (9) 

In Tables Ull and HVI we list for each of the diagrams the Mis needed for their evaluation. Masters denoted by an 
asterisk are of one-loop type. 




SE312m SE312md SE310m 




SE313m SE513m SE311m 



FIG. 5: The six two-loop 2-point Mis. External solid (dashed) lines describe on (off) -shell momenta. 



C. Two- loop 3-point Master Integrals 



There are nineteen two- loop 3-point Mis for Bhabha scattering. They are shown in Figures El to |H1 The equations 
for vertex Mis may be obtained with a differential operator: 



d 



1 



(2-s)p^ 



d 



(10) 



Another operator can be obtained by the exchange pi <-> P2- 

In Table IIVI we list all the vertex Mis for the various diagrams of Figure |21 The Mis for the evaluation of the 
QED (and QCD) vertex diagrams are worked out in [3ll [3^. The four Mis of prototypes V512ml and V512m2 are 
additionally needed for the evaluation of box diagrams Bl to B4 and have been determined in It was possible 
to determine all the vertex Mis with the method of differential equations. When dotted masters (or those with 
irreducible numerators) are involved one has to treat a system of coupled first order linear equations. For prototype 
V411ml, three masters had to be treated together. 



V614ml V614mld V614in2 

FIG. 6: The three two-loop vertex Mis with six internal lines. 




V512m2d V513m V513ind 

FIG. 7: The six two-loop vertex Mis with five internal lines. 




V414m V414md 



FIG. 8: 



The ten two-loop vertex Mis with four internal lines. 
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As an instructive example for the use of coupled differential equations we will study here the prototype V512m2. 
For few selected examples, we applied also alternative methods of calculation in order to have some cross checks; see 
also Two of them will also be exemplified here. One method (for MI V411m2) implies the integration over an UV 
divergent subloop and a subsequent subtraction in order to isolate an additional UV singularity, and the other one 
(for MI V512m2) introduces a subtracted dispersion relation. 



1. The master integrals for the prototype V512m2 



The prototype 

VR19 9^ h H h - -'—I [ d^fcld^fc2 (Plfcz)-'' {P2k2)-' 

has the masters V512m2 = V512m2(l,l,l,l,l;0,0) and V512m2d = V512m2(l,l,l,l,2;0,0), together with simpler ones. 
They contribute to the box diagrams Bl (planar double box) and B3 (nonplanar double box), see Table Hvl The 
momenta are, with the conventions used here, in the t channel. Applying the corresponding differential operator H42|l . 
in the t channel, we wrote the differential equations for the MI V512m2 and for V512m2N = V512m2(l,l,l,l,l;-1,0). 
One gets: 

t-^V512m2 = ■tiii_LM('v411m2[l,l,l,2;0l-V512m2 + 2 V512m2d^ 

dt 2{-4 + s + t)\ I ' ' ' ' J ; 

^ V512m2[l, 1,1, 1,2; 0,-1]- V512m2[l, 1, 1, 1, 2; -1, 0], (12) 



(-4 + s + <) 

t:^V512m2N = -tiii-LM ('2 V411m2fl, 1, 1, 2; -11 - (-2 + i) V411m2fl, 1, 1, 2; 01 

dt 4(-4 + s + i)V L > : > > J V y l > . > : j 

+ 4 V512m2N] + ^ V512m2[l, 1, 1, 1, 2; -1, -1] - V512m2[l, 1, 1, 1, 2; -2, 0]. (13) 

Additional functions at the right hand side are of prototype V411m2: 

-^^""^ r d^kid'^k2 (pife)"^ 

Hps + h+ /j2)2-m2]«[pi -p3-fci - k2]''[klY[t 



V411m2(a,6,c,d;<?) = --^ / , ,. , ' i,r^2icr^2i, ■ (14) 



Again, we write the basic master V411m2 (no dots, no numerators) without arguments. 

Besides Mis, at the right hand side of H12|) and H13|l one meets additional scalar integrals, which have to be expressed 
by algebraic equations through Mis: 

(l-5e + 6e2) / 2 - 3e \ 

V411m2 1,1,1,2;0 = -^^ ^ 4 V411m2 H SESllm , (15) 

^ ' ' ^ 2et \ -1 + 4e / ' ^ ' 

2-3e, , -l + 3e 
V411m2[l,l,l,2;-1] = (-2 SE310m + t SE311m) H V411m2, (16) 

(-l + e)2(-2 + 4e(l + e) + i(l-3e)) , e(l + 2e - t) 

V512m2 l,l,l,l,2;-2,0 = --^ ^ — ^ To i^Tlllm^ + ^V512m2 

^ ' ' ' ' ' ' ^ 8 l-2e 2e ^ -l + 2e 



+ [(8 - 2t - t^) + e(-56 + 18t + 6t^) + e^{U8 - 54i - llt^) + e^{-96 + 60t + 2t'^)] 

Im 

(-2 + 3e)(-2 + e(8 - 4t) + e^-A + t) + t) 



(-2 + 3e -2e + 4e2-t (-2 + t 

X , „ ^/ -T SESllm--^ — '-^^ ^V512ni2d 

16(-l + 2e -l + 4e 4e(-l + 2e) 



Set 

6t + t^ - 2e(-4 + <)(-5 + 2t) + 6^(48 + -20< + t^) 



8t 

Three further expressions may be left out here |0, but we have to give that for V512m2N 

4 , , t(l + 2e) 

V512in2N = (e V512m2) H ) ^V512m2d 

l-2e^ ' 2e(-l + 2e) 



SE310m 

V411m2 (17) 
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-SE311m+-!^ \ ^Tlllm^ 



8e(-l + 2e)(-l + 4e) 4(1 - 2e)2e 

2 + 9(-l + e)e -1 + e 

^SESlOmH V411m2. (18) 

AeH 4e ^ ' 

This equation may be inverted in order to eliminate V512m2d at the right hand side of (|12(l and 1)13(1 in favor of 
V512m2N, V512m2, and simpler Mis: 

V512m2N = -2 (e V512m2) - ^Tlllm^ + — — —7— r SESllm - SE310m 

dt ^ ' 2(-l + 2e) 2(-l + 2e)(-4 + t) ^ ' 

(6 - 16e - 2t + Set) 

- ^ 1^7—. ^ ^V411m2, (19) 

2(-4 + i) ' ^ ' 

d 2(l + 2e)-i -l + 2e,„ 
t— V512m2 = V512m2 H (2 V512m2N + V411m2) 

dt -4 + < -4 + < ^ ' 

+ 77^. TTT^. -SESllm-- ^ -Tlllm^. (20) 

2(-l + 2e)(-4 + <) (-l + 2e)(-4 + <) ^ ' 

The dependence on the unknowns in these two equations is such that they decouple as power series in e, and their 
solution is simplified thereby. Applying the inverted p8|) . the two masters are determined pol |: 

2x 

V512m2 = — (8C4 + 2C2i?[0, 0, x] - 4C2i?[0, 1, x] + H[0, 0, 0, 0, x] + 2H[0, 0, 0, 1, x]) + 0{e), (21) 

(—1 + a;^) 



V512m2d 



x x{l- H[0,x]-2H[l,x]) 



4e2(l-x)2 2e(l-x) 



2 



X 



■ [-4 - 15C2 - 4a; + C2a; + 4(1 + x) {H[0, x] + 2H[1, x]) 



4:{1-X)^{l+X) 

~ 2{x + 3) (iJ[0, 0, x] + 2iJ[0, l,x])- 8(1 + a;) (iJ[l, 0, x] + 2H[l, 1, a;])] + ©(e). (22) 

Another interesting opportunity appears when a system has some undotted and dotted masters being both UV 
finite, but IR divergent, and at the same time a related integral with numerator is both UV and IR finite [s^ . as 
it is the case here. Then, one may determine singularities in e from the algebraic relations between the functions; 
here from 1(18(1 . This way one also gets the divergent parts of 1(22(1 . This may serve as a cross check here. The dot in 
V512m2d introduces overlapping IR singularities and produces the l/e^ term. 

2. Example: The master integral V411m2 
The Ml V411m2 (see Figure (SI has a massless UV divergent subloop: 

V411m2 - ] [A:2][(fc, + fc2-pi)2][fc2-l][(A;i+p2)2]- ^^^^ 

We will use a subtraction procedure in order to isolate the remaining UV singularity. The two momentum integrations 
may be performed subsequently: 

d°k2 . c/2r(l-e)2r(e) 1 



[fc2][(fe + fci-pi)2] r(2-2e) [(A:i-pi)2]^' ^ > 

I - ^,, e(l + e)r(2e) 



idiv - / /r^-, ..\2 _ ^„^]2\2e- y"^^) 



dxdy x-^+'{l- x)^-'^' 
/o {[{1 - x){l - - xys]^) 

The Feynman parameter integral /div has a singularity at a: = and may be regulated by a subtraction: 

r-l /.I 



/div = / dxx-'+^{l-xY-^^ I d2/{[/(a;,y)-2^-/(0,2/)-2'] +/(0,y)-2^} 
JO Jo 



9 



r(£)r(2 - 2e) 
(l_4e)r(2-e) 



dx (1 - x) [x{l - xf 



\ f{x,y)-^^-f{0,y)-^^ 
ay : 



with 



fix, = (1 - x){l - yf - xys. 



(27) 
(28) 

(29) 



The remaining integrations in I^eg are regular and can be performed analyticaUy or numerically after the e-expansion: 



dx{\ — X 



5;]ln"-'=-V(x,y)lnV(0,2/) 



Lfe=o 



The first terms of the series expansion in e for V411m2 are (see Q): 

1 



V411m2 = 



2e2 2 



2(1 + 2;)^-' 2(1 + a;) 



(30) 



(31) 



They coincide with our results for the MI V411m2 given in |13l l2nl|. They are related, e.g., to the same MI given in 
[m by the following relation: V411m2 ^ 16[e<^''^r(l + (i)f[F-'^/{-2ef + F-'^/{-2e) + + . . ]. 



3. Example: The master integral V413md 



The master integral V413md is shown in Figure |S1 It is defined as follows: 

J — m2][{fc2 — (fci +pi)}2 — — m2]2[(fci — P2) 



V413md = — I , ^ m2 iriTTTS I3m7I TTrT' (32) 



The MI is UV and IR divergent, and one may first integrate over the UV divergent subloop and then treat the IR 
singularity by a subtraction. The first step gives |85j: 

V413md = ^/ [fc2„^2pp\„p^)2]5o [(pl + fer)^;m^m2]. (33) 

Due to the dotted photon propagator, the integral over fci produces for ki — > p2 an infrared singularity. This 
singularity may be isolated by a subtraction: 

V413md = So(.;™,m)B„(m;™,0) + -^y \kl - ^-^{{k, - ' 

(34) 

Here, the on-shell i?o function with underlined argument has a dotted line. The remaining integral over k\ is finite 
and one may evaluate it by use of a dispersion relation for the subtracted Bi^: 



1 /"^ 



1 1 



a-{pi+ki)^ a 



(35) 



with lmBQ{s,m^,m?) — tt-^/I — Am'^/s. After interchanging the order of integrations, one gets in the argument of 
the (7-integration two dotted, infrared-divergent one loop functions Co and Bq. In sum, the calculation of a two- loop 
vertex master integral has been reduced to the determination of one loop functions plus dispersion integrals: 



/•CO 

V413md = Bo {m^;m^,0) Bq (s;m'^,m^) - / — Imi?o(cr, m^, 



2 2. 2 n \ , ^0 {m^;rn^,0) 



Co {rn? , s, rr?\rr? , 0, cr) -f 



(36) 
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The dotted functions have been solved with the aid of IdSolver: 

D-2AQ{m'^) 



Bo (m2;m2,0) 
Co (to^, s, m^; m^, 0, tr) 



2 

e — — 7 Co (m^, s, TO^; m^, 0, fj) 



(37) 



rri^{(T — s) 



(1-6) 

(l-2e) 



ct(s " 2m2) - 2771^5 Ao(ct) 
m-^((7 — 4)(cr — s)^ cr 

:Bo(s;a,0) 



Ao(m2) 



2to2((T — 4to^)((t — s) TO^ 



m^((7 — s)' 



(ct - ^m?){a - s)) 



(38) 



The equations H37|l - H38|l are fulfiUed for any D, and the singularities of H36|l in e come only from the product of Bq 
functions. We show the first terms of the e-expansion for rri^ = 1: 



V413md = 



1 

2^ 



C2 



1+X 

1-x 



lii{x) 













1-x 



-Li2(-x) + ln(a;) ( 1 - ln(l - x) + ^ \n{x] 



da 



In ^'r^ - i V^W^) [ln(a) + 2 In (l - ^ 
2 2 L \ a 



0{e) 



(39) 



The result agrees with V413md given in the file MastersBhabha.m [Tjl, where it was derived with a differential equation. 
The first two coefficients are easily read from (|39|l . and the constant term is: 



V413ind[0, x] 



(2 



l-x 4(1 -a;) 



1 _l_ 77* 1 4- T 

ln^(a;) + —— {\n{x) [1 - ln(l - x) - 31n(l 



1-x 



- 3Li2(-a;) - Li2(a;)} , 



(40) 



D. Two- loop 4-point Master Integrals 

There are thirty three two-loop box Mis. They are shown in Figures Elton] Several analytical expressions for box 
Mis in terms of HPLs may be found in the file MastersBhabha.m [13. In Table im we list the double-box diagrams 
to which they contribute. Additionally, one may see the other contributing masters. By now, few of the two-loop 
box Mis are known analytically. For diagram Bl, the Mis B714ml and B714mlN are given in [l3|; for diagram B2, 
the MI B714m2 is known as two-dimensional integral representation [l^; for diagram B3, the leading divergent part 
of MI B714m3 is published The technique used is the Mellin-Barnes representation [s^. Issl l56l| in combination 
with summation techniques a la j57l. 58, 59] . Other box Mis were derived with systems of differential equations. For 
diagrams Bl and B3 the MI B512ml has been given in j23|. For diagram B5, a set of two double-box Mis, B514m 
and B514mN, was derived recently with the restriction to electrons and photons (one flavor) in F'or our set of 

masters, we use instead of B514mN a dotted function: 



B514md 



-xy 



(l-x2)(l-2/2) 



H[0, x]H[0, y] ^ _ _ ^^^^^^ ^ 



e 

- 4C2i?[0, y] + G[-l/y, 0, x]H[0, y] + G[-y, 0, x]H[0, y] - 2H[0, y]H[-l, 0, x] 
-6H[0, x]H[-l, 0, y] + 4iJ[0, x]H[0, 0, y] + G[-l/y, x] [3(2 + H[0, 0, y]] 

-G[-y, x] [3C2 + [0, 0, y]] + 2H[0, x]H[l, 0, y] - H[0, 0, 0, y]] +0{e). (41) 

The expression depends on HPLs H, but also on two-dimensional HPLs G (see Appendix A of [g^I; we use here the 
notations of and the G[—y, 0, 0, x] as given in l20|). 

A differential operator for the derivation of differential equations for box diagrams in the s-channel is: 



d If,. s{P2-P3)} d 



The corresponding operator in the t-channel is obtained by the replacements s t and p2 <-> — Ps- Since the number 
of differential operators p'^d/dp'^ is larger than the number of kinematic invariants, there is some freedom of choice. 
Let us note that representation l|42|) is much simpler than e.g. that chosen in 113 • 
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B714in3 B7141m3dl B714m3d2 

FIG. 9: The nine two-loop box Mis with seven internal lines. 



E. Dotted Mis and Mis with irreducible numerators 



From the point of view of automatic calculations, there is no essential difference whether a set of Mis is given with 
or without Mis with irreducible numerators. There exist algebraic relations to transform from one set to the other. 
However, when the Mis are determined with the help of differential equations, the solution may come out much easier 
when the unknown Mis are properly chosen. Introducing numerators or dots on lines will change the nature of a 
MI concerning both the UV and IR singularities. Related to this is that some coupled differential equations for Mis 
can be decoupled due to separated orders of their e-expansions. Moreover, the solution for a given MI can be much 
simpler compared to that for another choice. As an example may serve (14 I J, which is by far more compact than the 
corresponding MI B514mN with numerator, given in Equations (37)-(39) of 1 181. Let us also note that the singularity 
of B514mN is absent in l|4H) . Contrary, the masters of prototype SE312m [lj| are an example for the opposite case 
where a solution with numerator is simpler than that with a dotted propagator. 

Let us finally mention that any of the dotted Mis can be replaced by a MI with an appropriate numerator. However, 
it is not always possible to "move" a dot arbitrarily from one line to another line. E.g. the dot in V512m2d cannot be 
moved from the massless line to one of the two internal massive lines. The arising integral V512m2 [2, 1, 1, 1, 1;0, 0] 
(defined in Hll|l 'l may be considered as a master integral, but cannot replace V512m2d. This may be easily proven by 
the following two algebraic relations: 

V512m2 1, 1, 1,2, 1:0.0 V512m2d ^ z -^V411m2 

^ ' ' ' ' ' 2x(l + 2e) 

3(-l + 2e)(-2 + 3e)(-l + 3e) 3(2 - 13e + 27e2 - ISe^) , ^ 

+ — 7^^^^ ^SE311m+-^ ^ ^SESlOm, (43) 

42;(-l + 2e + 8e2) ^ 2ex^{l + 2e) ' ^ ' 

l-5e + 6e2 (2 - 13e + 27e2 - ISe^) 

V512m2 2,1,1,1,1;0,0 = -e V512m2 V411m2 + -^^ ^ ^SESllm. (44) 

2ex 4e.i;(l — 4e) 

We conclude for this specific example that the set of masters could be chosen to contain one master out of the pair 
of integrals with a dot on a massless line ( V512m2d, V512m2 [1, 1, 1, 2, 1;0, 0]), and one master out of the pair of 
integrals with no dot or a dot on a massive line ( V512m2, V512m2 [2, 1, 1, 1, 1;G, 0]). 



1. The divergent parts of the master integrals for the prototypes B512m2 and B512m3 



As was explained in Section IIII C II under certain conditions one may determine singular parts of Mis in a purely 
algebraic way or with a combination of algebraic relations and differential equations. We have used this method in 



B613ml 



B613mld 




B613m2d 




B613m3d2 





B613m3 




B613m3d3 



B613m2 




B613m3dl 




B613m3d4 



B613m3d5 

FIG. 10: The ten two-loop box Mis with six internal lines. 



order to determine the singularities of the so far unknown MI for prototypes B512in2 and B512m3: 

1 X 



B512m2 ^ 



B512m3d2 



e x^ — \ 



B512m2dl = - 



B512m2d2 = 
B512m3 = 

B512m3dl = 



(l-a;)2 

£2 8a:(l - yf 
1 V 



(4C2 + H[Q, 0, x] + 2H[Q, 1, x]) + 0(1), 



^ + -{2 + H[fd,x\ + 2H[l,x\) +0{l), 



-(l + a;2)(l + y)3 + C2[8(l-a;)2(l-y)y] 



e 4a;(l - y)2(l + yf 
+ 2(1 + y) {x(\ - yf + 2y (1 + x^)] (if [0, y] + 2iJ[l, y]) 

- (1 - x^){l + yfH% x] + 2j;(l - xf{l - y) (i7 [0, 0, y] + 2ff [0, 1, y])] + 0(1), 
1 xyi7[0, a;] 

"^(l-.T2)(l-y)2 

1 2a;y 
" 7(l-x2)(l-t/)2 



iJ[l, 0, x] - [-1, 0, x] + iJ[0, 0, x] 



I + [0, x] (i? [0, 2/] + 2 iJ[l, y])]+ 0{l), 



1 

B512m3d3 = 



e2(l-y2) 



[4C2 + i?[0,0,j/] + 2i/[0,l,y]] + O(l). 
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B512m3d2 B512m3d3 



B513m 



B513mdl 



B513md2 




B513md3 B514m B514md 

FIG. 11: The fourteen two-loop box Mis with five internal lines. 

F. Additional master integrals with several flavors 



In QED, there are additional ferniion flavors besides electrons. In two-loop Bhabha scattering this leads to the 
additional diagrams SE5f, V4f, B5f shown in Figure^] They may be derived from Figures to O by a replacement of 
the closed electron loop in diagrams SE5, V4, B5 by a loop with a second mass scale M. As a consequence, additional 
Mis will arise. We have determined them with DiaGen/IdSolver, and Table IVl lists the Mis which contribute to the 
evaluation of Feynman integrals of the new prototypes. 



IV. SUMMARY 



The results presented here document a further step towards a complete two-loop prediction of small angle Bhabha 
scattering, needed e.g. for a precise luminosity determination at the future International Linear Collider. 

We have determined a complete set of scalar master integrals needed for the calculation of the virtual two-loop 
corrections to massive Bhabha scattering in QED, both for one flavor and for several flavors. The Mis are shown 
pictorially, and we tabulate to which prototype diagrams they will contribute. Further, we determined the flrst terms 
of the e-expansions for all the two-loop self-energy and vertex master integrals for the one-flavor case, and some of the 
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SE312Mln 



SE312Mlmd 



V412M2m 






V412M2md 



V412Mlm 



V412Mlmd 



B512M2md B512M2m 
FIG. 12: The eight additional master integrals with two different mass scales. 



two-loop box masters. Some of these master integrals were unknown so far. The analytical expressions are collected 
in the Mathematica file MastersBhabha.m and are publicly available As by-products of our numerous tests of 
the results presented here, we expressed HPLs up to weight 4 by a minimal basis (Mathematica file HPL4.m ^3)' 
have generalized the sector decomposition algorithm for the evaluation of Feynman integrals in the Euclidean region 
to the case of irreducible numerators. 

Once the last, most complicated master integrals are determined one will have to combine the virtual two- loop 
corrections to Bhabha scattering, together with their counter terms, with electroweak corrections and with a package 
for the treatment of real bremsstrahlung. 
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TABLE I: The number of (two-loop + one-loop) master integrals needed to calculate the two-loop vertex diagrams and box 
diagrams (Fig. 



Diagram 


VI 


V2 


V3 


V4 


V5 


Bl 


B2 


B3 


B4 


B5 


B6 


tadpole MI 


0+1 


0+1 


0+1 


0+1 


0+1 


0+1 


0+1 


0+1 


0+1 


0+1 


0+1 


2-point Mis 


3+1 


4+1 


4+1 


1+1 


3+1 


4+1 


5+2 


5+1 


4+2 


3+2 


3+2 


3-point Mis 


4+0 


10+0 


5+0 


2+0 


1+0 


7+0 


11+1 


13+0 


10+1 


4+1 


4+1 


Box type Mis 












9+0 


15+0 


22+0 


11+1 


2+1 


3+1 


total 


7+2 


14+2 


9+2 


3+2 


4+2 


20+2 


31+4 


40+2 


25+5 


9+5 


10+5 


net 


9 


16 


11 


5 


6 


22 


35 


42 


30 


14 


15 



15 



TABLE 11: The one-loop and two-loop Mis needed for the evaluation of the self-energy diagrams 



MI 












rcf. 


Tlllm* 


+ 


+ 


+ 


+ 


+ 


r /I /1 1 

I44| 


SE212m* 


+ 


+ 








m 


SESllm 






+ 






oms:[47l 


SE312m 


+ 


+ 








[32, 52] 


SE312md 


+ 


+ 








[52j,sec. inrni 


SE313m 








+ 


+ 


oms:f47. 4^. 491 



TABLE III: The 1-point, 2-point and 3-point Mis entering basic two- loop vertex and box diagrams in Figures to |3 Stars 
denote one-loop Mis. 



MI 


Bl 


B2 


B3 


B4 


B5 


B6 


VI 


V2 


V3 


V4 


V5 


rcf. 


i iiim 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 




SE210m* 




+ 




+ 


+ 


+ 












IMi 


SE212m* 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


[A Al 

[44J 


SESlOm 


+ 




+ 


















13^>^>^ 










+ 




1 

+ 


+ 


1 

+ 


+ 




+ 


oms. |4 ( ] 


SE312m 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 




+ 


Too KOl 












1 


1 




1 








[^^91 ^an ITTT Til 

^zj,oec. 1111 rj 


SE313m 




+ 


+ 


+ 


+ 






+ 


+ 


+ 




oms: [47, 48, 49] 


SE513m 




+ 






















V311m* 




+ 




+ 


+ 


+ 












[44] 


V411ml 




+ 




+ 




+ 














V411ml[dl— d2] 




+ 




+ 




+ 














V411m2 


+ 




+ 


















[31J 


V412ml 


+ 




+ 








+ 


+ 










V412m2 


+ 


+ 


+ 


+ 




+ 


+ 


+ 


+ 




+ 


m 


V413m 




+ 


+ 


+ 


+ 






+ 


+ 






[32] 


V413md 




+ 


+ 


+ 


+ 






+ 


+ 






[32] 


V414m 




+ 


+ 


+ 


+ 






+ 


+ 


+ 




m 


V414md 




+ 


+ 


+ 


+ 






+ 


+ 


+ 




[32] 


V512ml 




+ 




+ 
















[20] 


V512mld 




+ 




+ 
















[20] 


V512m2 


+ 




+ 


















m 


V512m2d 


+ 




+ 


















m 


V513m 


+ 




+ 








+ 


+ 








[32] 


V513md 


+ 




+ 








+ 


+ 








m 


V614ml 






+ 










+ 








[32] 


V614mld 






+ 










+ 










V614m2 




+ 






















total = 25-1-4* 


11+2* 


16+4* 


16+4* 


14+4* 


7+4* 


7+4* 


7+2* 


14+2* 


9+2* 


3+2* 


4+2* 
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TABLE IV: 4-point Mis entering basic two-loop box diagrams in Figure|3 An asterisk denotes one-loop MI. Mis with a dagger 
are not known completely analytically; see text. 



MI 


Bl 


B2 


B3 


B4 


B5 


B6 


ref. 


R71 Ami 














[171 


R71 Am 1 M 
















R71 Am9 
















R71 AmO Phi H*^! 




1 












R71 Am^ 






1 

-r 








n Qit 


rS ( I'HIlo LCli — uzj 






1 










B613ml 
















Rfil ^m1 H 


1 




1 










Rfil ^mO 




1 




1 








B613m2d 




+ 




+ 








Rfil ^m^ 
DO ±01110 
















Rfil ^m^ Phi 






r 

-r 










PCI Om-I 


1 

+ 




1 

+ 








[901 






+ 








1 

+ 


oec. liii rj 11 


pen 1 ^oT 








+ 






oec. liii rj ii 


pen Om^ 
D0X.^II10 


1 

-r 




1 












-L 

r 




I 








Sec HIT Elf 


B513m 




+ 


+ 


+ 








B513m [dl— d3] 




-1- 


+ 


+ 








B514m 




+ 


+ 


+ 


+ 




[12] 


B514md 




+ 


+ 


+ 


+ 




Sec. IlIIDl 


B412m* 








+ 


+ 


+ 




total = 33-f 1* 


9 


15 


22 


u+r 


2+r 


3+r 





TABLE V: The two-loop master integrals for several flavors. An asterisk denotes a one-loop MI. 



MI 


SE5f 


V4f 


B5f 


ref. 


Tlllm* 


+ 


+ 


+ 


M 


SE212m* 




+ 


+ 


[44] 


SE210m* 






+ 


[44J 


SESllm 






+ 


oms: f47l 


SE312m 






+ 


[32, 52] 


SE312md 






+ 


[^,Sec. HTlBl 


SE312Mlm 


-1- 


+ 


+ 


oms:f6l1 


SE312Mlmd 


+ 


+ 


+ 


oms: [611 


V311m* 






+ 


[44] 


V412Mlm 






+ 




V412Mlmd 






+ 




V412M2m 




+ 


+ 




V412M2md 




+ 


+ 




B412m* 






+ 


[1^,44] 


B512M2m 






+ 




B512M2md 






+ 
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APPENDIX A: NUMERICAL EVALUATION OF MULTI-LOOP DIAGRAMS WITH IRREDUCIBLE 

NUMERATORS IN THE EUCLIDEAN REGION 

Here we derive expressions for the Feynman integrals G{X) defined in Q with L loops, a numerator function X{ki), 
N internal lines, and E external legs with momenta (see also e.g. 62]). The propagator momenta qi compose as: 

= c\ki + ... c'^kL + ... + dipi + ... + cI'ePe- (A1) 

In [Hill, a method was derived for the numerical evaluation of multi-loop Feynman integrals in the Euclidean region. 
The method is based on sector decomposition (see also ^3, 64]). In [TslllGl ]. the method was worked out explicitely 
for propagators only, including dotted ones (i^i ^ 1), i.e. for X = 1. We need the formulae also with numerators and 
repeat the starting expression for completeness here, in our notations: 

r(x) ^ ^ f d^k,...d^k^ X 

The denominator of G{1) contains, after introduction of Feynman parameters Xi, the momentum dependent function 

N 

^2 - ^a;,(gf-m2) = kMk-2kQ + J. (A3) 



m 

1=1 

The linear terms in may be eliminated by a shift: 

r-i 



k = k + M-'Q, (A4) 
and the matrix M in ljA3|l may be diagonalized by a rotation: 

k^k' = V k, (A5) 

kMk = k'Md^agk' (A6) 

- J^a.fcf, (AT) 

Md^ag = MV-^ = (ai,...,aL). (A8) 

After the diagonalization, the momentum integrals may be easily performed, and the remaining Feynman parameter 
integral contains the term: 

(det A/)-^/2 _ jj^^)N^-D{L+l)/2 

(^2)^--^i/2 " ^(a.)A'„-DL/2 ' ( 9) 

with 

A*' = -(J-QM-^Q), (AlO) 
= vi + ...VN, (All) 

and the definitions [stI ]: 

U{x) = (det M), (A12) 
F{x) = (det M)^? = - (det M) J + QMQ. (A13) 

The M is defined with — M /det M. Finally, one arrives at: 

Equation (3) in contains formula (jA14|) for i^i — 1. 

From the above it is straightforward to formulate the parameter Feynman parameter integrals for the case of 
nontrivial numerators. The corresponding formulae for the simplest cases are: 



U{x) 



N^-l-D{L + l)/2 



Fix 



N^-DL/2 



J2 ^~IuQ 



I 

- a 



(A15) 
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and 



G{ki^k2l3) = (-1) 



E 



T{vi) . . .V{vn) Jo fJi 



1 N 




[MllQlUM2lQlh 



U{x)F{x) 



Fix)- 



ai 



(A16) 



The general case is treated in |65j- Potential UV singularities arise from the overall factors with F functions and 
from U{x), while IR singularities arise from F{x). The Mathematica code sectors. m performs the isolation of 
IR singularities by sector decomposition and the integrations in the Euclidean region. For the numerical calculations 
the package CUBA [g^I is used. 



APPENDIX B: HARMONIC POLYLOGARITHMS UP TO WEIGHT 4 



The master integrals in the file Master sBhabha.m jl^ are expressed by harmonic polylogarithms (HPLs) up 
to weight 4. Harmonic polylogarithms fulfill simplified algebraic relations of harmonic sums. One may use this and 
determine a basis of functions by direct evaluation, and then express the others by algebraic relations. For details we 
refer to and references therein. 

The three HPLs of weight 1 are: 



H[0,x] = ln(a;), 
H[l,x] = -ln(l-x), 
H[-l,x] = ln(l + a;). 



(Bl) 
(B2) 
(B3) 



The nine HPLs of weight 2 may be expressed by three independent HPLs of weight 2 plus those of lower weight: 



H[0,l,x] - Li2(a;), 
H[0,^l,x] ^ -Li2(-a;) 

H\-\,\,x\ = Li2 ' 



Lis 



In 2 H{--\,x\. 



(B4) 
(B5) 

(B6) 



The other HPLs are determined with the relation i/[a, 6, x\ — —H[b, a, x] + H[a, x]H[b, x]. 

There are 27 HPLs of weight 3 [H3> and 8 of them have to be added to the basis: one for each of the six index sets 
of the class {a,a,b) (see equations (3.22), (3.23) of 



H[0,0,l,x] 
H[0,0,-l,x] 
F[0,l,l,x] 
H[0,^l,-l,x] 

H[-l,l,l,x] 
ff[l,-l,-l,x] 



Li3(a;), 
-Li3(-x), 
'S'i,2(a;), 
Si.2{-x), 

^..^ 

. 1-x 

-i>i,2 



- Si,2 



In 2 
+ ln2 



Li2 
Li2 



1 + a; 
2 

1-x 



Lis. 



iln2 2i7[-l,x], 



1 



ln^2 H[l,x] 



and two for the class (a, 6,c) (see equations (3.7) to (3.10) of p^): 

1 



H[l,0,-l,x] 



5i,2(a;^) - H[0, 1, 1, x] - [0, -1, -1, x] + H[l, x]H[0, -1, a;] 



H[Q,-i,i,x] = 

Jo V 



Lio 



l + y 



Li2 ( I 



In 2 H[0,-l,x\. 



(B7) 
(B8) 
(B9) 
(BIO) 

(Bll) 
(B12) 



(B13) 

(B14) 



Further, H\a , a, a, x] = ^H[a, x]"^. The integral in ljB14p may be performed analytically, see V3{x) in (|B37|) or equation 
(188) in [ia. 
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There are 81 HPLs of weight 4: the three index sets of class {aaaa) have 1 element each (none is basic), the six 
index sets of class {aaah) have 4 elements each (one is basic), the three index sets of class (aabb) have 6 elements each 
(one is basic), and the three index sets of class {aabc) have 12 elements each (three are basic). The 18 independent 
HPLs of weight 4 may be chosen to be: 

i/[0, 0,0,1,2;] = huix), (B15) 
iJ[0, 0,0,-1,2;] = -Li4(-a;), (B16) 
ff[0,l,l,l,x] = ^1,3(2;), (B17) 

H[-l, l,l,l,x] = -Li4 (^) + Li4 + ln(l - x)Us (^] - hn'il - x)U, ^ ^ ~ ^ 



^ 2 J \2J ' ' \ 2 J 2 

- iln='(l-x)lni±^, (B18) 

i7[0,-l,-l,-l,a;] = -5i,3(-a;), (B19) 

H[l,-l,-l,-l,x] = if[-l,l,l,l,-x], (B20) 

H[0,0,l,l,x] = 52,2(x), (B21) 

i?[0,0,-l,-l,x] = 52,2(-a;), (B22) 



i7[l,l,-l,-l,x] = / J^H[l,-l,-l,y 

J a 1 - ?y 



= -A'i(a;) +ln2 J2(a;) + iln^2 i/[l,l,a;], (B23) 

i/[-l,0,0,l,a;] = / J^Usiy), (B24) 
Jo 1 + 2/ 

H[-l, 0, 1, 0,x] ^ - ^5*2,2 (x^) - H[-1,Q, 0, 1, x] + H[0, 0,-1, -1,2;] + H[0, 0, 1, 1, a;] 
- H[~l,x]H[0,0,l,x] 

+ H[0,x]{-H[l,x]H[0,~l,x] + H[-l,x]H[0,l,x] + H[l,0,-l,x]), (B25) 



ij[-i,i,o,o,x] = - r-^[5i,2(i-y)- ^1,2(1)] 

Jo 1 + y 



1 r dy 2r 



= - ln(l + x) [5i,2(l - x) - 5i,2(l)] - - / In^(y) ln(l + y), (B26) 

2 Jo 1 - y 

iy[0,-l,l,l,2;] = C/4(a;)-ln2 iJ[0,-l,l,a;]- iln2 2 i/[0,-l,a;], (B27) 

iJ[0, 1,-1,1,2;] = -F4(2;)-2iJ[0,-l,l,l,2;] + ln2 {H[1,0, -1, x] - H[l,x]H[0, ~l,x]) , (B28) 

iJ[0, 1,1,-1,2;] = 14(2;) + ^W^4(a;) +iJ[0,-l,l,l,2;] -ln2 (i7 [1, 0, -1, 2;] - iJ[l, 2;]iJ[0, -1, 2;]) , (B29) 

i?[0, 1,-1, -1,2;] = -J74(-.x) + ln2 i/[0, 1,-1,2;] - iln2 2 i7[0,l,2;], (B30) 

^[0,-1,1,-1,2;] = \4(-a;)-2iJ[0, 1,-1,-1, x]- In 2 (i7 [1, 0, -1, 2;] - iJ[l, 2;]iJ[0, -1, 2;]) , (B31) 

i?[0, -1,-1, 1,2;] = -Vi{-x)-^Wi{~x) + H[0,l,~l,-l,x] 

+ In 2 (i7[l,0,-l,2;] -i?[l,2;]i?[0,-l,2;]). (B32) 

Further, H[a,a,a,a,x] = ■^H[a,x]'^. Some auxiliary integrals are defined below in this appendix. The integrals Ki 
and J2 in IIB23p may be performed analytically, see ljB39p and (|B40p , and that in (|B26p is not known to us in analytical 
form. We just mention that the functions with index sets (0000), (0001), (0011), (Gill), (1111) may also be found 
in [b^I • Our list of HPLs is available in file HPL4 . m in ^] . We give here one explicit example of a HPL of the class 
(aabc), derived with equation (4.26) of [H [Hj: 

i7[-l,l,0, 1,2;] = 2i/[0, 1,1,-1, 2;] + iJ[0, 1,-1,1,2;] +i7[l,.T](-iJ[l, -1,0,. t] - iJ[l,0,-l,x] 

+ i7[-l,0,l,2;] +iJ[0,-l,l,x]) - 2i7[-l,2;]iJ[0, 1,1,2;] +i/[l,-l, x]i/[l, 0,2;]. (B33) 

The above list of functions is a quite compact basis which is sufficient for the computation of all HPLs until weight 4. 
For comparisons with [sij we used additionally to HPLs with arguments x also the following functions: 

Linix^) = 2"-i[iJ[0, 0,1,2;] -iJ[0,..., 0,-1, x]] (B34) 
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with (n-1) zeroes in the arguments, and: 

Siaix^) = H [0, 1, 1, x^] = 2 [i? [0, 1, 1, x] - H[0, -1, 1, x] - H[0, 1, -1, x] + H[0, -1, -1, x]] , (B35) 
S2,2{x^) = H [0, 0, 1, 1, x^] = 4 [ff [0, 0, 1, 1, x] - H[0, 0,-1, 1, a;] - H[0, 0, 1,-1, x] + iJ[0, 0, -1,-1, x]] . (B36) 

For a study of boundary conditions, one needs the values of HPLs at a; = and x = 1. Here we follow the algorithms 
proposed in '28"|. 

We have introduced several auxihary functions. For H[0, —1, l,x] (see (|B14|) ') we need: 



dy 
y 



Ji(x) +ln(a;) 



1 + 2/ 



Li2 



- Li2 ( i 



1+ X 



Li2 - 



In 2 [Li2(-a;) + ln(x) ln(l + x)] , 



(B37) 



with 



- dy 



= ^Lis 

+ ln 
1 



1 + y 

1 - X 



h\{y) ln(l - y) 



2 

2x 

l-x 



- Lis {-x) - Lis - 



l-x 



Li3 - 



Li2 



1 - X 



Li2 



2x 

l-x 



-ln(l-x) 



- In 2 ln(l - x) + 2 In 2 ln(2a;) + 2 ln(a;) In 



2x 

l-x 

+ In (2a:) Li2 (-a;) + In 
1 + a; 



The i?[l, 1, —1, — l,a:] (see (jB23|l ') is known analytically, using: 



1 - X 



Li2 



1 - a; 



(B38) 



Ki{x) 



dy 

1-y 



i-y 



In 2 J3(a;)--i^2(x)-ln(l-.T) 



dy 
1 - y 



Li2 



Jz{x)-\n{l-x) 



Li2 



Lio - 



1 - a; 



Li2 



l-x 



-^1.2 1 2 



1 



-\t?2\t?(1- x), 

4 ^ ' 



■ln21n'^(l - a;) 



(B39) 



(B40) 



with 



■h{x) 



K2{x) 



dy 



= -Lis 



1 - y 

1 - .a: 



ln(l - y) ln(l + y) 



Li, ( \ 



ln(l — a:)Li2 



1-a: 



dy 
1 - y 



7 



In(l-y) ln2(l + y) 



1 



-iln21n2(l-a;), 



1 



= — C2 - C2 In^ 2 + In'* 2 + - ln^(l + x) - ln-'(l + a;) ln(l - x) + - ln^(l + x) ln^(l - a:) 

fl — a;) 

- In 2 ln(l - a;) ln(l + a;) In - In^ 2 In (l - a;^) + (In^ 2 + ln^(l - a;))Li2 [(1 - a;) 12] 

2(1 + x) 

+ ln2((l - x)/{l + x))Li2[-(l - x)/{l + x)] + (In' 2 - W{1 + x))Li2[(l + x)/2] 

-2 ln(l - x)Li3[(l - a;)/2] - 2 ln[(l - x)/{l + a;)]Li3[-(l - a;)/(l + a;)] 

+2 ln(l + a:)Li3[(l + a;)/2] + 2Li4[(l - a;)/2] + 2Li4[-(l - a;)/(l + a;)] - 2Li4[(l + x)/2\. 



(B41) 



(B42) 



In the class (0,1,1,-1) we have an analytical expression for (iJ[0, 1, 1, — 1, a;] + i?[0, 1, — 1, 1, a;]), but not for 
i?[0, — 1, 1, 1, a;] and (i?[0, 1, 1,-1, x] — if[0, 1,-1, 1, a;]). An analogous statement applies for the class (0,-1,-1, 1). 
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The reason is that W4{x) is known from equations (A. 28) to (A. 33) of reference |51| : 



V 

1 



ln"(l-y) ln(l + y) 



15 



Si,^{x') + 2Si,3{-x) - [Li4 {lo) - Li4 {-uj)] + — C4 

+ \n{Lo) [Lia {uj) - Lig (-tj)] - \ \v?{uj) [Lis (w) - Liz (-w)] - i \n^{Lo) \n{x), 

2 b 



(B43) 



with u! = (1 — x)/(l + x). By expressing Si^^ix'^) through HPLs with argument x 28], one may derive the remarkably 
simple relation 



Wiix) + Wii-x) = ~Si^3{x^) + 2Si,3ix) + 2Si,3i-x). 



(B44) 



This relation follows also from |5l|. 

We have rewritten the two remaining auxiliary functions: 



V 



Si, 



1 + 2/ 



ln2 Ji(a;)--/3(x)+ln(a;) Su 

1 o 

In 2 [ln(x) ln(l + x) + Li2(-x)] 



1 + a; 



Si.:. 



(B45) 



Viix) 



^ dy 



ln(l - v) 



y 

hix) - Li2{x) 
1 



Li2 
Li2 



1 + y 
2 

l + x 



- Li2 

— Li2 



ln2 



Si,2{x') - Si^2{x) - 5i,2(-a;) + ln(l + x)Li2(x) 



We mention finally that we do not know analytical expressions for the following auxiliary functions: 

dy 



h{x) 

h{x) 
h{x) 



1 + y 

" dy 

1 + y 

dy 

1 + y 



Li3(y), 
Li2(y) ln(l-2/), 



dy ,„2 



1 + 2/ 



\n\y)\n{l~y). 



(B46) 

(B47) 
(B48) 
(B49) 
(B50) 



The first integral is related to 1, 0, 0, 1, x] (see (IB24ll ). and the next two integrals to four combinations of the six 
HPLs of classes (0, a, a, 6); see (jB27|l to (jB32|) and the remark related to before (|B43|) . The last one is related to 
H[-l,l,0,0,a;], see (IB^ . 
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